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Abstract 

The  effect  of  gas  diffusion  medium  (GDM)  intrusion  on  the  performance  of  proton  exchange  membrane  (PEM)  fuel  cells  is  investigated.  The 
mechanical  behaviors  of  various  GDM  are  characterized  in  compressive,  flexural,  and  shear  tests.  The  results  are  used  in  a  numerical  model  to 
calculate  the  channel  intrusion  of  GDM.  The  intrusion  calculation  from  the  numerical  model  agrees  well  with  the  measurements  from  an  intrusion 
measurement  setup  and  a  pressure  drop  measurement  device  for  various  GDM.  A  simplified  reactant  flow  redistribution  model  of  parallel  channels 
developed  in  this  study  suggests  that  a  5%  variation  in  GDM  intrusion  can  result  in  a  20%  reduction  of  reactant  flow  in  the  most  intruded  channel. 
The  GDM  intrusion  and  intrusion  variation  are  found  to  induce  significant  performance  discrepancy  among  cells  of  a  30-cell  automotive  fuel  cell 
stack  consisting  of  two  different  production  lots  of  commercial  GDM.  The  study  suggests  that  in  the  future  mass  production  of  fuel  cell  stacks,  GDM 
manufacturers  need  to  greatly  tighten  their  product  variations  in  mechanical  property  and  thickness  to  ensure  reliable  PEM  fuel  cell  operations. 
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1.  Introduction 

A  proton  exchange  membrane  (PEM)  fuel  cell  stack  is  com¬ 
prised  of  a  series  arrangement  of  repeating  cell  units,  each  of 
which  consists  of  a  membrane  electrode  assembly  (MEA)  where 
the  electrochemical  reaction  takes  place;  two  layers  of  gas  diffu¬ 
sion  media  (GDM)  for  the  distribution  of  fuel  and  oxidant  gases 
over  the  catalytic  surfaces;  and  a  bipolar  plate  which  directs  fuel 
and  oxidant  gases  from  the  inlet  manifolds  through  a  network  of 
flow  channels.  During  the  fuel  cell  stack  assembly,  these  repeat¬ 
ing  cells  are  stacked  together  with  intermediate  seal  gaskets  and 
then  compressed  to  provide  adequate  gas  sealing,  as  well  as  to 
reduce  contact  resistances  at  the  material  interfaces.  The  GDM 
plays  an  important  role  in  PEM  fuel  cells  by  (1)  acting  as  a 
buffer  layer  between  the  softer  MEA  and  stiffer  bipolar  plates; 

(2)  acting  as  diffuser  for  reactant  gases  traveling  to  electrodes; 

(3)  transporting  product  water  to  the  gas  flow  channels;  (4)  con¬ 
ducting  electrons;  (5)  transferring  heat  generated  at  MEA  to  the 
coolant  within  the  bipolar  plates.  GDM  is  normally  composed 
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of  porous  non- woven  carbon  fiber  paper  or  carbon  fiber  cloth. 
The  effect  of  GDM  compression  on  the  fuel  cell  performance  has 
been  investigated  by  many  researchers  mostly  on  electrical,  ther¬ 
mal,  and  mass  transport  resistances  within  the  bulk  material  and 
across  the  contact  interfaces  [1-9].  To  minimize  contact  resis¬ 
tance,  high  contact  pressure  between  the  lands  of  reactant  gas 
flow  field  and  gas  diffusion  media  is  usually  desired.  However, 
these  studies  generally  found  that  optimal  compression  pres¬ 
sure  would  need  to  be  determined  by  balancing  its  conflicting 
effect  on  GDM’s  porosity,  diffusivity,  permeability,  electrical 
conductivity,  and  thermal  conductivity. 

It  is  worth  noting  that  most  papers  on  the  fuel  cell  compres¬ 
sion  and  its  effect  on  the  PEM  fuel  cell  performance  considered 
only  a  relatively  small  domain  for  the  purpose  of  understanding 
material  behavior.  These  included  a  half-channel-and-half-land 
domain  in  computational  modeling;  ex  situ  clamping  setups  for 
samples  of  several  cm2  in  size;  and/or  small-scale  single  cells 
of  5  to  50  cm2.  When  considering  the  compression  in  full-scale 
(stacks  with  100  or  more  cells)  PEM  fuel  cell  applications,  the 
flow  distribution  factor  starts  playing  a  significant  if  not  dominat¬ 
ing  role  [10].  Flow  distribution  in  channels  of  manifold  devices 
is  a  topic  inherited  from  plate  heat  exchangers  comprising  dis¬ 
crete  minichannels,  for  example,  with  inlet  and  outlet  manifolds 
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Nomenclature 

a 

width  of  GDM  in  the  shear  test 

b 

thickness  of  GDM  under  compression  in  the  shear 
test 

A  b 

displacement  of  shear  piston  in  the  shear  test 

bo 

initial  thickness  of  GDM  in  the  shear  test 

d 

channel  depth 

D 

hydraulic  diameter  of  flow  channel 

GDM 

gas  diffusion  media 

/ 

GDM  channel  intrusion 

MEA 

membrane  electrode  assembly 

n 

number  of  parallel  channels 

P 

force  applied  on  the  shear  piston  in  the  shear  test 

PEM 

proton  exchange  membrane 

Q 

flow  rate  per  channel 

t 

time 

f^avg 

average  velocity  of  reactant  gas  flow 

W 

channel  width 

Greek  symbols 

<5 

normalized  Q  decrease 

£ 

normalized  D  decrease 

P 

reactant  viscosity 

opening  on  one  side  (U-configuration)  or  on  opposite  sides  (Z- 
configuration)  [11,12].  In  recent  years  several  papers  on  flow 
distribution  in  fuel  cells  and  fuel  cell  stacks  were  published.  Bar- 
rears  et  al.  [13]  reported  on  flow  velocity  distribution  measured 
in  a  single  50  cm2  cell.  Grega  et  al.  [14]  studied  velocity  distri¬ 
bution  in  a  scaled  up  model  of  a  21 -cell  stack.  Maharudrayya  et 
al.  [15,16]  modeled  flow  distribution  for  a  variety  of  single  cell 
flow  field  configurations.  Koh  et  al.  [17]  proposed  a  systematic 
algorithm  for  computation  of  pressure  and  flow  distribution  in 
internal  gas  manifolds  of  a  100-cell  stack.  Because  the  number 
of  cells  is  relatively  large,  the  flow  in  manifolds  has  to  be  treated 
as  turbulent  whereas  in  most  previous  publications  both  channel 
and  manifold  (header)  flows  are  considered  laminar.  In  ah  these 
papers,  deviation  from  uniform  channel-to-channel  or  cell-to- 
cell  flow  distribution  is  attributed  mainly  to  the  geometry.  The 
results  may  serve  as  a  design  guide.  For  example,  Park  and  Li 
[18]  considered  the  effect  of  flow  and  temperature  distribution 
on  the  performance  of  a  5 1 -cell  stack,  and  one  of  the  conclusions 
is  that  the  hydraulic  diameters  of  the  inlet  and  outlet  manifolds 
should  be  at  least  20  times  larger  than  the  channel  diameter. 

A  more  comprehensive  approach  to  flow  distribution  was 
implemented  by  Chang  et  al.  [19].  A  conventional  finite  element 
approach  to  flow  branching  and  merging  back  along  the  inlet  and 
outlet  headers  was  coupled  with  reactant  consumption  and  water 
production  along  the  channels,  which  defined  the  actual  cell  per¬ 
formance.  The  model  was  validated  against  experimental  data 
for  a  190-cell  Ballard  Mk  7  stack,  and  the  fitted  friction  factors 
were  averaged  over  the  length  of  the  channel  and  header.  The 
most  distinctive  part  of  this  paper  is  the  case  study  of  so  called 
deviation  from  manufacturing  tolerances,  which  was  achieved 


by  introducing  friction  coefficient  variation  for  a  single  cell,  or 
by  imposing  some  distribution  of  friction  coefficients  on  ah  100 
cells  of  the  model  stack.  By  introducing  this  simulated  pertur¬ 
bation,  the  authors  proved  that  a  relatively  small  deviation  may 
cause  significant  changes  of  flow  in  a  single  cell  or  the  entire 
cell  to  cell  flow  distribution  in  the  stack. 

Due  to  the  compression  difference  in  the  GDM  over  chan¬ 
nels  and  lands,  the  unsupported  portion  of  the  GDM  can  intrude 
into  reactant  gas  channels  upon  assembling  fuel  cell  stacks.  The 
direct  consequence  of  GDM  intrusion  is  that  it  can  increase  the 
pressure  drop  of  the  reactant  gases  in  the  intruded  channels. 
Nitta  et  al.  [7]  reported  the  measurement  of  GDM  intrusion  at 
various  degrees  of  compression  and  channel  widths  using  a  dial 
indicator  in  an  intrusion  measurement  setup.  They  reported  that 
the  intrusion  increased  as  the  degree  of  compression  increased. 
However,  the  sensitivity  to  channel  width  was  much  less  signifi¬ 
cant,  which  could  be  caused  by  very  thick  and  soft  GDM  used  for 
testing.  In  a  subsequent  paper  by  Hottinen  et  al.  [8]  the  effect  of 
inhomogeneous  land  versus  channel  compression  was  taken  into 
account  in  a  fuel  cell  model,  which  showed  significant  impact 
on  the  fuel  cell  performance.  In  this  study,  because  the  model 
domain  was  limited  to  single  channel/channel,  the  difference  in 
performance  was  mostly  attributed  to  the  under-the-land  condi¬ 
tions  (contact  resistance  and  permeability),  while  the  possibility 
of  channel-to-channel  reactant  flow  maldistribution  due  to  GDM 
intrusion  was  not  considered  at  ah.  More  recently,  Basu  et  al. 
[20]  investigated  the  two-phase  reactant  flow  maldistribution 
resulted  from  GDM  intrusion  caused  by  uneven  compression 
distribution  over  the  active  area  of  a  single  cell.  The  authors  also 
gave  a  more  comprehensive  explanation  for  the  implications  of 
reactant  flow  maldistribution.  It  was  concluded  that  poor  reac¬ 
tant  distribution  not  only  can  cause  performance  and  stability 
problems,  but  also  create  conditions  for  catalyst  corrosion,  and 
therefore  affect  durability. 

Until  recently,  relatively  little  attention  has  been  paid  to 
studying  GDM’s  mechanical  behavior  and  its  effect  on  fuel  cell 
performance.  Mathias  et  al.  [1]  reported  the  mechanical  charac¬ 
terization  of  GDM.  The  compressive  behavior  was  characterized 
by  placing  a  GDM  between  two  flat  plates  and  the  deflection 
was  measured  as  a  function  of  compressive  force.  It  was  found 
that  after  the  first  loading,  the  GDM  generally  showed  differ¬ 
ent  unloading  and  reloading  curves,  suggesting  weakening  of 
the  material  from  the  first  loading.  To  characterize  the  bending 
behavior  of  GDM,  Mathias  et  al.  proposed  the  use  of  flexural 
tests  such  as  ASTM  D790.  With  the  load/deflection  response 
from  the  flexural  test,  flexural  modulus  and  flexural  strength 
could  be  determined.  Lai  et  al.  [21]  investigated  the  compressive 
pressure  between  GDM  and  ME  A  over  the  lands  and  channels 
in  a  single-channel  cell  and  in  fuel  cell  stack.  They  found  that 
fuel  cells  could  lose  significant  compression  pressure  due  to 
the  hysteresis  in  GDM’s  compressive  behavior.  Because  of  the 
relatively  small  ratio  between  the  channel  width  and  the  GDM 
thickness,  typically  ranging  from  2  to  10,  the  shear  behavior 
of  GDM  can  also  affect  the  compressive  pressure  distribution. 
From  finite  element  modeling  and  testing  using  PressureX™ 
(Sensor  Product  Inc.)  film,  they  suggested  that  a  higher  trans¬ 
verse  shear  modulus  was  preferred  in  the  GDM  to  increase  the 
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contact  pressure  over  the  channels  and  to  reduce  the  potential 
ME  A  buckling. 

In  this  paper,  an  effort  to  develop  a  GDM  intrusion  model 
based  on  the  finite  element  method  (FEM);  experiments  to  mea¬ 
sure  the  GDM  intrusion  and  flow  field  pressure  drop  to  validate 
the  intrusion  model;  and  a  simplified  reactant  flow  redistribu¬ 
tion  model  to  estimate  the  effect  of  channel  intrusion  on  the 
reactant  flow  redistribution  will  be  reported.  How  PEM  fuel  cell 
performance  can  be  affected  by  GDM  channel  intrusion  and 
intrusion  variation  as  a  direct  consequence  of  GDM’s  property 
and  thickness  variations  will  also  be  discussed. 

2.  GDM  channel  intrusion  model 

In  this  study,  a  2D  finite  element  model  was  used  to  deter¬ 
mine  the  intrusion  of  GDM  into  the  gas  channel.  The  model  was 
built  using  the  commercial  finite  element  software,  ABAQUS® . 
Fig.  1  illustrates  the  deformed  shape  of  the  model,  which  consists 
of  a  single  channel  flow  field  and  a  single  GDM  layer  con¬ 
strained  at  the  lower  surface.  The  GDM  intrusion  is  defined  as 
the  vertical  distance  from  the  apex  of  GDM  surface  facing  the 
reactant  gas  channel  to  the  plane  of  lands.  The  land  and  channel 
profile  is  modeled  as  rigid  surface.  Contact  condition  is  consid¬ 
ered  between  the  rigid  surface  and  the  GDM  upper  surface.  Two 
element  types,  a  gasket  and  a  plane  strain  element,  are  super¬ 
posed  together  and  used  to  model  the  GDM.  Three  mechanical 
tests,  a  compressive  test,  a  flexural  test,  and  a  shear  test,  were 
conducted  to  determine  the  material  properties  to  be  used  in 
the  model.  It  should  be  noted  that  there  is  yet  no  standardiza¬ 
tion  of  methods  to  measure  the  GDM’s  mechanical  properties. 
Although  a  more  comprehensive  analytical  and  experimental 
effort  would  be  needed  to  completely  characterize  the  GDM 
mechanical  properties,  it  is  found  that  the  tests  used  in  this  study 
can  satisfy  our  needs  to  calculation  the  GDM  intrusion. 

Following  Mathias  et  al.  [1],  the  GDM  compressive  behavior 
was  characterized  by  placing  a  GDM  between  two  flat  circular 
steel  blocks.  The  load  and  deflection  relationship  was  measured 
to  obtain  the  compressive  stress-strain  curve  in  the  thickness 
direction.  The  in-plane  elastic  behavior  was  characterized  by  a 
flexural  test  described  in  ASTM  D790  from  which  the  modu¬ 
lus  of  elasticity  was  determined.  To  measure  the  shear  property 
of  GDM,  a  punch-and-die  type  of  shear  test  was  developed. 
In  the  test  fixture  schematically  shown  in  Fig.  2(A),  a  GDM 
sample  was  clamped  between  two  sets  of  rectangular  clamping 
blocks  separated  by  a  fixed  distance.  A  rectangular  shear  pis¬ 
ton  was  pushed  down  on  the  unsupported  GDM  to  generate  a 
load-defection  (P  vs.  A b)  curve.  An  effective  transverse  shear 
modulus  was  then  determined  by  using  the  initial  linear  portion 


Fig.  1.  The  schematic  of  GDM  intrusion  of  a  GDM  compressed  by  a  single 
channel  flow  field. 


(B) 


Fig.  2.  Shear  test  setup  and  stress-strain  schematic  for  the  transverse  shear 
modulus  of  GDM. 

of  the  load-deflection  curve  as  seen  in  Fig.  2(B).  It  should  be 
noted  that  the  flexural  and  shear  tests  used  in  this  study  may  pro¬ 
vide  only  an  approximation  to  the  Young’s  modulus  and  shear 
modulus  of  the  GDM.  However,  as  will  be  seen  in  the  later  dis¬ 
cussion,  the  good  agreement  between  the  intrusion  measurement 
and  the  model  calculation  for  several  different  GDM  types  sug¬ 
gests  that  these  tests  can  generate  useful  material  information 
for  model  calculations. 

In  the  finite  element  model,  the  compressive  stress/strain 
curve  from  the  compressive  test  was  used  to  model  the  behav¬ 
ior  of  the  gasket  elements.  The  solid  elements  were  assigned 
orthotropic  elastic  properties  with  a  very  small  through-plane 
modulus.  Moduli  of  elasticity  from  the  flexural  test  were  used 
as  Young’s  modulus  in  the  in-plane  directions.  Apparent  shear 
moduli  from  the  shear  test  were  used  as  the  shear  moduli.  Since 
GDM  has  a  very  high  open  porosity  (commonly  at  80%),  a  very 
small  Poisson’s  ratio  was  used. 

To  understand  the  GDM  intrusion  behavior,  a  carbon  paper- 
based  GDM  denoted  as  “GDM  X”  was  used  as  the  model  GDM. 
The  thickness  of  the  GDM  is  about  0.260  mm.  The  effective 
transverse  shear  modulus  from  the  shear  test  is  about  14.2  MPa 
in  the  first  in-plane  direction  (denoted  as  “Dir.  1”)  and  9.2  MPa 
in  the  second  in-plane  direction  (denoted  as  “Dir.  2”).  The  mod¬ 
uli  of  elasticity  determined  from  the  flexural  test  are  900  MPa 
and  300  MPa  in  Dir.  1  and  Dir.  2,  respectively.  The  measured 
compressive  strains  in  the  thickness  directions  are  tabulated  in 
Table  1. 

Fig.  3  illustrates  the  intrusion  versus  compressive  pressure 
over  the  land  for  a  sheet  of  GDM  X  compressed  by  a  flow  field 
of  1  mm- wide  channel  and  2  mm- wide  land  against  a  rigid  plane. 
Note  that  the  GDM  intrusion  in  this  configuration  is  equivalent 
to  the  one  in  a  single  fuel  cell  with  a  symmetrical  plane  at  the 
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Table  1 


Compressive  stress  and  strain  data  of  GDM  X 


Strain  (mm  mm  1 ) 

Stress  (MPa) 

0.15 

0.7 

0.21 

1.4 

0.27 

2.1 

0.33 

2.8 

0.35 

3.1 

0.37 

3.45 

mid-plane  of  the  MEA.  In  Fig.  3,  the  line  denoted  “GDM  X — Dir 
1”  represents  the  case  where  the  first  in-plane  direction  corre¬ 
sponds  to  the  channel  width  direction.  At  a  compressive  pressure 
of  2.8  MPa  over  the  land,  the  intrusion  is  0.073  mm  in  Dir.  1  and 
0.079  mm  in  Dir.  2.  For  the  whole  compression  range,  it  is  seen 
that  the  intrusion  for  the  case  with  lower  stiffness  in  both  shear 
and  bending  (Dir.  2)  is  higher  than  the  case  with  higher  stiff¬ 
ness  (Dir.  1),  which  illustrates  that  stiff er  GDM  can  have  less 
intrusion. 


Fig.  5  illustrates  schematically  a  typical  set  of  data  collected 
using  the  intrusion  test  fixture.  The  difference  between  X  and  Y 
under  the  compression  load  of  interest  gives  the  absolute  GDM 
intrusion  into  the  slot.  A  relative  intrusion  can  also  be  determined 
by  (X  —  Y)/X. 

To  validate  the  model  calculation,  a  sheet  of  GDM  X  was 
compressed  at  2.8  MPa  over  the  land  and  the  intrusion  was 
measured.  Fig.  6  illustrates  the  relative  intrusion  comparison 
between  the  model  calculation  and  experimental  results.  Good 
agreement  is  seen  for  both  in-plane  directions. 

Fig.  7  illustrates  similar  comparisons  for  five  types  of  GDM 
from  four  different  manufacturers.  Good  agreement  is  also  seen 
for  all  cases  in  spite  of  a  various  degree  of  compressibility,  bend¬ 
ing  stiffness,  and  shear  moduli,  which  is  reflected  by  a  relatively 
wide  spread  of  the  relative  intrusion  ranging  from  about  75%  to 
more  than  95%.  The  good  agreement  between  the  experiment 
and  model  calculation  further  suggests  that  the  GDM’s  intrusion 
behavior  can  be  reasonably  captured  using  the  mechanical  tests 
and  the  intrusion  model  described  in  this  paper. 


3.  GDM  channel  intrusion  experiment  and  model 

validation  4.  Reactant  flow  redistribution  model 


To  validate  the  intrusion  model  calculation,  an  intrusion  mea¬ 
surement  setup,  as  schematically  shown  in  Fig.  4,  was  developed. 
In  this  intrusion  fixture,  a  slot  coupon  was  machined  to  mimic 
the  geometry  of  single  channel  flow  field.  A  floating  bar  was 
placed  in  contact  with  the  GDM  in  the  slot.  Two  extensometers 
were  mounted  at  the  two  ends  of  the  floating  bar  to  monitor 
the  movement  of  the  floating  bar  when  compression  load  was 
applied  from  the  top  of  the  fixture.  Two  tests  were  performed. 
The  first  test  was  carried  out  by  placing  a  rigid  shim  on  top  of  the 
GDM  as  shown  in  Fig.  4(A).  The  extensometers  were  zeroed  at 
35  kPa  when  the  load  was  applied  from  the  top  of  the  fixture.  As 
the  compression  load  was  increased,  the  displacement  of  floating 
bar  was  recorded  as  X ,  which  represents  the  thickness  change 
of  the  GDM  over  the  lands.  A  second  test  as  shown  in  Fig.  4(B) 
was  performed  in  a  similar  way  by  monitoring  the  GDM  thick¬ 
ness  change  in  the  slot  under  compression.  In  the  second  test, 
the  rigid  shim  was  not  used,  and  the  movement  of  the  floating 
bar  was  recorded  as  Y,  which  was  also  zeroed  at  35  kPa. 


In  order  to  estimate  the  effect  of  GDM  intrusion  into  flow 
field  channels  on  flow  distribution,  an  n  parallel-channel  con¬ 
figuration  is  considered  as  shown  in  Fig.  8,  which  consists  a 
set  of  n  parallel  channels  of  the  same  length  with  a  nominal 
hydraulic  diameter  D  and  a  nominal  flow  rate  per  channel  Q. 
Using  a  simplified  model  based  on  Hagen-Poiseuille  equation 
for  incompressible  laminar  flow  in  cylindrical  conduits  [23]: 


dp  _  ^/x£/ayg 
dv  D2 


(1) 


where  p  is  the  pressure  of  reactant  gas;  v  is  the  distance  from  the 
inlet;  /z  is  the  viscosity  of  the  reactant  gas;  £/av g  is  the  average 
velocity,  all  the  channels  would  have  the  same  pressure  drop  A p 
from  inlets  to  outlets,  given  by 


(2) 


APPLIED  GDM  COMPRESSION  PRESSURE  OVER  LANDS  (MPa) 

Fig.  3.  GDM  intrusion  vs.  the  compressive  pressure  over  the  lands  of  a  model 
GDM  X  for  the  case  of  1  mm- wide  channel. 


To  investigate  the  flow  redistribution,  the  case  where  a  chan¬ 
nel  has  more  GDM  intrusion  than  the  nominal  channels  is 
considered.  Due  to  this  additional  intrusion  the  hydraulic  diam¬ 
eter  of  the  more  intruded  channel  is  decreased  by  AD  compared 
to  the  rest  of  the  channels.  In  this  case  the  flow  through  this 
channel  will  decrease  by  an  amount  of  A Q.  Because  the  system 
is  maintaining  a  constant  flow  rate,  the  amount  A Q  will  be  re¬ 
distributed  among  the  remaining  (n  —  1)  channels.  Provided  that 
the  pressure  drops  in  the  nominal  and  the  more  restricted  chan¬ 
nels  are  still  equal  and  the  change  in  channel  hydraulic  diameter 
of  the  more  intruded  channel  is  relatively  small,  the  following 
equation  can  be  derived: 

Q  +  AQ/(n  —  1)  Q-AQ 

D4  (D  -  AD)4  ^  ’ 
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Fig.  4.  Schematic  of  the  GDM  intrusion  test  setup. 


Fig.  5.  Schematic  of  load-displacement  curves  obtained  from  intrusion  test. 


Dividing  Eq.  (3)  by  Q/D 4,  and  substituting  8  =  AQ/Q,  and 
s  =  AD/D,  Eq.  (3)  can  be  rewritten  in  a  dimensionless  form  of 


1  + 


8 

(n  -  1) 


0-3) 
0  -  £)4 


(4) 


Solving  Eq.  (4)  relative  to  8  =  (1  —  8)  the  flow  in  the  more 
restricted  channel  as  a  percentage  to  the  nominally  expected 
flow  is  given  by 


8(n,  s)  = 


n(  1  —  £)4 
(n  —  1)  +  (1  —  e)4 


(5) 


Dir.1  Dir.2 

Fig.  6.  Comparison  in  relative  intrusion  between  the  model  calculation  and 
experimental  measurement  in  GDM  X  samples  of  two  different  orientations. 


Eq.  (5)  is  presented  in  Fig.  9  for  the  cases  of  2,  5,  15,  and  100 
channels.  As  one  would  expect,  in  the  case  of  100  channels, 
the  excessive  flow  coming  from  the  single  restricted  channel  is 
much  more  easily  distributed  between  the  remaining  99  channels 
than  the  case  of  2  channels.  Therefore,  the  flow  rate  in  the  more 
restricted  channel  changes  with  the  magnitude  of  the  restriction 
much  more  significantly  for  the  case  with  more  channels. 

Note  that  this  simplified  flow  redistribution  model  is  just  to 
provide  an  estimate  for  a  header  fed  set  of  parallel  channels, 
which  is  aimed  at  demonstrating  the  strong  correlation  between 
GDM  intrusion  and  flow  distribution.  Similar  to  Chen  et  al. 
[22],  where  the  flow  distribution  in  the  PEM  fuel  cell  stack  was 
estimated  without  accounting  for  any  electrochemistry  or  heat 
transfer,  we  consider  isothermal  dry  reactant  distribution.  In  this 
case,  if  assuming  uniform  compression  over  the  active  area  and 
perfect  reactant  gas  distribution  in  the  manifolds,  the  approach 
can  be  easily  applied  to  fuel  cell  stacks  in  which  GDM  intrusion 
from  channel  to  channel  within  each  cell  is  uniform  but  not  from 
cell  to  cell,  i.e.  some  of  the  cells  will  have  higher  than  nominal 
hydraulic  resistance.  For  example,  this  can  happen  in  fuel  cell 
stacks  where  a  small  batch  of  thicker  or  softer  GDM  is  mixed 
with  a  larger  batch  of  thinner  or  stiffer  GDM,  where  intrusion  is 
similar  within  the  same  batch  of  GDM  but  significantly  different 
between  the  two  batches.  According  to  the  100-channel  case  in 
Fig.  9,  a  stack  comprising  hundreds  of  cells  that  are  operated 
at  relatively  high  reactant  utilization  conditions  (i.e.  at  a  low 
stoichiometry),  a  5-15%  additional  decrease  in  hydraulic  diam¬ 
eter  can  cause  a  20-50%  decrease  in  flow  in  the  cells  with  the 
most  severe  GDM  intrusion.  Therefore,  if  the  stack  is  operated 
at  relatively  high  reactant  utilization  conditions  (i.e.  at  low  sto¬ 
ichiometry),  then  the  cells  with  the  highest  GDM  intrusion  will 
experience  reactant  starvation,  which  can  lead  to  failure. 

It  should  be  once  again  noted  that  the  flow  redistribution 
model  described  in  this  paper  represents  a  greatly  simplified 
model.  Many  aspects  of  this  model  can  be  further  refined  to  bet¬ 
ter  reflect  the  reality  of  a  rather  complicated  problem.  Ultimately, 
a  two-phase  flow  3D  CFD  model  accounting  for  electrochemi¬ 
cal  reaction  should  be  a  better  choice  if  required.  However,  the 
simplified  model  presented  in  this  paper  is  believed  to  serve  the 
purpose  for  an  initial  understanding  of  how  GDM  intrusion  vari¬ 
ation  may  contribute  to  the  problems  in  fuel  cell  performance.  It 
is  also  believed  to  provide  useful  information  for  the  direction  of 
further  improvement  in  materials  and  flow  field  designs.  As  will 
be  seen  in  the  later  discussion  on  the  performance  of  a  30-cell 
fuel  cell  stack,  the  model  is  useful  enough  to  explain  observed 
cell-to-cell  performance  variations. 
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GDM  SAMPLES 

Fig.  7.  Relative  intrusion  comparison  between  experimental  and  model  calculation  for  a  wide  range  of  GDM  types  and  manufacturers. 


Fig.  8.  Schematic  of  a  header  fed  parallel-channel  configuration. 


5.  Results  and  discussions 

Since  channel  intrusion  reduces  the  hydraulic  diameter  of 
reactant  gas  channels,  a  direct  effect  of  channel  intrusion  is 
the  increase  of  the  reactant  gas  pressure  drop  from  the  inlet 


to  outlet  of  a  flow  field.  Thus,  pressure  drop  measurement 
can  also  provide  an  estimate  of  channel  intrusion  and  intru¬ 
sion  variation,  and  vice-versa.  To  measure  the  pressure  drop, 
a  series  of  flow  fields  with  various  land  and  channel  designs 
were  fabricated  and  tested.  The  data  to  be  discussed  in  Fig.  10 
were  from  a  flow  field  that  was  composed  of  6  straight  chan¬ 
nels  of  300  mm  long,  0.4  mm  wide,  and  0.26  mm  deep  with 
a  land  width  of  2  mm.  In  each  measurement,  a  single  piece 
of  GDM  was  clamped  between  this  flow  field  and  a  flat  plate 
with  an  average  of  1.46  MPa  compression.  Air  flow  rate  of 
1  slpm  was  introduced  at  the  inlet  and  pressure  drop  was  then 
measured.  To  investigate  the  effect  of  GDM  intrusion  varia¬ 
tion  on  pressure  drops,  a  number  of  samples  from  several  lots 
of  another  commercial  gas  diffusion  media,  dubbed  GDM  Y, 
from  several  production  lots  were  used.  These  samples  exhibit 
a  range  of  mechanical  property  variation  in  compressibility 
(compressive  strain  measured  at  1.2  MPa  of  compressive  pres¬ 
sure,  slightly  lower  than  1.46  MPa  used  in  the  pressure  drop 


Fig.  9.  Reactant  flow  in  the  cell  with  additional  blockage  compared  to  the  nominal  cells. 
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Fig.  10.  Comparison  of  pressure  drop  measurement  and  model  calculation  for  GDM  samples  representing  various  mechanical  property  variations. 


test),  shear  modulus,  and  modulus  of  elasticity.  Fig.  10  illus¬ 
trates  the  pressure  drop  versus  compressive  strains  of  various 
GDM  samples  measured  at  1 .2  MPa  of  compressive  pressure. 
The  compressive  strain  from  different  lots  of  GDM  at  1.2  MPa 
has  a  range  of  0.212-0.268.  It  is  seen  that  the  measured  pressure 
drop  increases  from  44  kPa  to  54  kPa  as  the  compressive  strain 
increases. 

Although  it  is  possible  to  directly  calculate  the  pressure  drop 
analytically,  e.g.,  from  a  CFD  program,  it  has  not  been  attempted 
since  the  current  interest  is  in  the  relative  change  in  pressure 
drop  for  GDM  with  different  intrusion.  Instead,  the  hydraulic 
diameters  of  the  intruded  channels  are  estimated  from  the  well- 
known  hydraulic  diameter  equation  of  a  rectangular  flow  channel 
[23] 


2  w{d  —  I) 
(w  +  d  —  I) 


(6) 


where  w  is  the  channel  width;  d  is  the  channel  depth;  I  is  the 
GDM  intrusion.  Noting  from  Eq.  (2)  that  the  pressure  drop  at 
a  fixed  flow  rate  is  inversely  proportional  to  the  fourth  power 
of  hydraulic  diameter,  one  can  use  the  experimentally  measured 
pressure  drop  for  the  GDM  of  compressive  strain  of  0.218  at 
1.2  MPa  as  a  baseline  and  then  calculate  the  pressure  drops  for 
other  GDM  samples  using  Eq.  (7): 

D4  (baseline  GDM) 

A  P(new  GDM)=  A  /^baseline  GDM)  — - -  (7) 

D4(new  GDM) 

It  is  seen  in  Fig.  10  that  the  calculated  pressure  drop  agrees 
very  well  with  the  measured  pressure  drop,  which  further 
suggests  the  usefulness  of  the  channel  intrusion  model  and  sim¬ 
plified  flow  model.  The  wide  variation  in  pressure  drop  observed 
for  various  lots  of  GDM  Y  also  demonstrates  that  pressure  drop 
is  indeed  very  sensitive  to  the  GDM  channel  intrusion  varia¬ 
tion,  which  is  the  consequence  of  GDM  property  and  thickness 
variations. 


Fig.  11.  Reactant  gas  flow  in  fraction  of  nominal  flow  in  the  more  intruded  channel  for  four  GDM  with  increasing  deviation  in  compressive  strain,  thickness,  shear 
modulus,  and  modulus  of  elasticity  at  various  channel  depths. 
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Fig.  11  illustrates  the  sensitivity  of  the  reactant  gas  flow, 
expressed  in  fraction  of  the  nominal  flow,  in  the  channels  with 
the  additional  GDM  intrusion  compared  to  the  nominal  chan¬ 
nels.  A  flow  field  of  1  mm- wide  channel  and  2  mm- wide  land 
is  used  for  this  study.  Four  cases  are  considered  with  increasing 
deviations  in  compressive  strain,  thickness,  shear  modulus,  and 
modulus  of  elasticity  from  Case  1  to  4,  which  are  tabulated  in 
the  insert  of  Fig.  11.  In  addition,  four  different  channel  depths 
of  0.25  mm,  0.5  mm,  0.75  mm,  and  1  mm  are  considered.  In  this 
example,  the  intrusion  is  calculated  using  the  intrusion  model 
by  applying  a  2.76  MPa  of  compressive  pressure  over  lands  with 
GDM  X.  The  intrusion  is  then  used  to  estimate  the  reduction  in 
hydraulic  diameter  from  Eq.  (6).  With  the  normalized  hydraulic 
diameter  change  s  determined,  the  normalized  flow  rate  in  the 
restricted  channel  is  calculated  using  Eq.  (5).  In  the  nominal 
case,  denoted  “Case  1”,  the  GDM  X  is  oriented  so  that  Dir.  1 
is  in  the  width  direction  of  lands  and  channels,  which  results  in 
a  GDM  intrusion  of  0.073  mm.  As  expected  in  Fig.  1 1  the  flow 
in  the  nominal  channels  is  100%  for  all  channel  depths.  In  Case 
2,  a  compressive  strain  deviation  of  +0.06  at  2.76  MPa  is  intro¬ 
duced  to  the  channel  of  interest.  At  a  channel  depth  of  1  mm, 
the  deviation  of  reactant  gas  flow  for  the  cell  using  this  GDM  is 
about  5.6%.  As  the  channel  depth  decreases,  the  flow  deviation 
increases.  At  a  channel  depth  of  0.25  mm,  the  flow  deviation 
is  about  27%.  In  Case  3,  a  thickness  deviation  of  +0.015  mm 
is  further  added  to  Case  2  for  the  channel  of  interest,  which 
shows  a  gas  flow  deviation  of  7.6%  at  1  mm-deep  and  35.2%  at 
0.25  mm-deep  channels.  In  Case  4,  a  shear  modulus  deviation 
of  —5.4  MPa  and  modulus  of  elasticity  deviation  of  —600  MPa 
are  added  to  Case  3.  Case  4  is  analogous  to  a  situation  such 
that  a  few  cells  in  a  100-cell  stack  use  Dir.  2  of  GDM  X  in  the 
channel  width  direction  while  the  rest  of  the  cells  use  Dir.  1 .  In 
this  case,  one  can  find  a  flow  deviation  of  10.5%  at  1  mm-deep 
and  46.3%  at  0.25  mm-deep  channels.  Note  that  the  deviation  in 
compressive  strain,  thickness,  and  shear  modulus  illustrated  in 
Fig.  11  are  not  unusual.  In  fact,  they  all  have  been  observed  in 
the  commercial  GDM  X  from  different  production  lots  in  our 
study. 

Figs.  10  and  11  provide  an  intriguing  glimpse  into  how 
GDM  properties  are  tightly  related  to  flow  distribution  and  flow 
field  geometry.  This  implies  that  GDM  manufacturers  have  to 
improve  the  variation  to  a  much  tighter  specification.  For  fuel 
cell  producers  in  general,  there  is  a  continual  drive  to  operate 
fuel  cells  at  as  low  stoichiometry  as  possible  in  order  to  increase 
the  efficiency  of  fuel  cell  systems.  Therefore,  the  threshold  of 
performance  instability  is  decreasing.  A  simple  example  is  that 
if  a  fuel  cell  stack  operates  at  1 . 1  stoich  and  the  flow  rate  in 
one  of  the  cells  drops  below  90%  of  the  nominal  flow,  the  per¬ 
formance  in  that  cell  would  become  unstable.  Thus,  to  use  the 
GDM  X  with  the  given  variation  found  in  Case  4,  the  flow 
field  would  have  to  be  designed  with  a  channel  depth  larger 
than  1  mm  or  the  fuel  cell  would  have  to  be  operated  at  a  sto¬ 
ich  much  higher  than  1.1.  Neither  option  is  desirable  when 
the  fuel  cell  industry  pursues  an  increasingly  higher  volumetric 
power  density  and  reactant  gas  utilization.  Therefore,  it  is  impor¬ 
tant  that  GDM  manufacturers  can  produce  GDM  of  superior 
mechanical  properties — higher  compressive,  shear,  and  bending 


stiffness  while  significantly  tightening  thickness  and  property 
variations. 

To  investigate  the  GDM  intrusion  effect  on  automotive  fuel 
cells,  a  30-cell  fuel  cell  stack  which  contained  equal  number 
of  GDM  Y  from  two  different  production  lots  was  tested.  For 
convenience,  the  two  lots  are  denoted  as  Lots  A  and  B  and  the 
corresponding  cells  are  called  Lot  A  and  Lot  B  cells.  In  the  sin¬ 
gle  fuel  cell  testing  with  50  cm2  active  area  of  graphite  flow  field 
and  a  channel  depth  of  about  1  mm,  both  lots  of  GDM  produced 
nominally  the  same  polarization  curves  under  the  same  fuel  cell 
operating  condition,  suggesting  that  they  had  the  same  thermal, 
electrical,  and  mass  transport  properties.  The  only  noticeable 
difference  between  the  two  production  lots  were  in  the  mechan¬ 
ical  properties  and  thickness.  Under  the  stack  compression,  the 
intrusion  of  Lot  A  GDM  is  determined  to  be  0.069  =b  0.006  mm. 
The  intrusion  of  Lot  B  is  0.078  d=  0.009  mm.  It  should  be  noted 
that  the  intrusion  and  its  variation  are  calculated  based  on  18 
samples  for  each  lot.  The  intrusion  value  suggests  that  the  Lot 
A  GDM  is  generally  stiffer  than  Lot  B.  Furthermore,  the  intru¬ 
sion  variation  in  Lot  A  is  tighter  than  Lot  B,  suggesting  that  the 
mechanical  properties  in  Lot  A  GDM  are  more  consistent  from 
sample  to  sample.  The  pressure  drop  measurements  as  discussed 
for  Fig.  10  also  support  intrusion  calculation  since  the  44kPa 
pressure  drop  in  Lot  A  GDM  is  much  lower  than  that  of  5 1  kPa 
in  Lot  B,  also  suggesting  that  Lot  A  GDM  has  less  intrusion 
than  Lot  B.  In  the  subsequent  fuel  cell  testing,  the  30-cell  stack 
was  tested  under  various  cathode  stoichiometry  at  a  current  den¬ 
sity  of  1.2  A  cm-2  to  investigate  the  cell  voltage  sensitivity  to 
cathode  stoichiometry. 

Fig.  12  illustrates  the  cell  voltage  of  Lot  A  and  B  cells  nor¬ 
malized  with  respect  to  the  averaged  cell  voltage  at  a  cathode 
stoichiometry  of  1.8.  It  is  seen  that  for  all  cathode  stoichiometry, 
Lot  A  cells  have  higher  mean  cell  voltages  than  those  of  Lot  B 
cells,  which  correlates  well  with  the  lower  calculated  intrusion 
and  pressure  drop  measurement.  As  the  cathode  stoichiometry 
is  decreased,  the  separation  of  mean  cell  voltages  between  Lots 
A  and  B  increases.  Although  the  cell  voltage  variation  in  both 
lots  of  cells  increases  as  the  cathode  stoichiometry  is  decreased, 
the  variation  increase  in  Lot  B  cells  is  much  more  pronounced 


Fig.  12.  Cell  voltage  variation  at  various  cathode  stoichiometry  for  two  different 
production  lots  of  GDM  Y. 
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than  in  Lot  A.  Finally,  as  the  cathode  stoichiometry  dropped 
below  1.3,  the  Lot  B  cells  became  so  unstable  that  the  stack 
ceased  to  operate.  To  determine  the  performance  anomaly  in 
a  fuel  cell  stack  is  often  a  difficult  task  as  the  number  of  the 
root  cause  is  often  more  than  one.  In  this  study,  it  is  believed 
that  the  GDM  intrusion  has  played  an  important  role  leading  to 
the  observed  stoichiometry  sensitivity  because  the  stack  perfor¬ 
mance  data  have  correlated  well  with  the  GDM  intrusion  data 
in  many  aspects  which  can  be  summarized  as  follows:  (1)  The 
higher  mean  cell  voltage  in  Lot  A  cells  than  Lot  B  could  be 
caused  by  the  lower  intrusion  and  hence,  the  higher  flow  rate 
in  Lot  A  cells.  (2)  Since  Lot  B  cells  have  worse  intrusion  and 
intrusion  variation,  the  flow  maldistribution  in  Lot  B  cells  would 
make  them  increasingly  sensitive  to  the  decreased  cathode  flow, 
which  leads  to  increasing  spread  in  cell  voltage  as  the  cathode 
stoich  is  decreased.  (3)  Finally,  in  the  few  most  intruded  Lot  B 
cells,  the  reduced  overall  cathode  flow  rate  could  bring  the  flow 
rate  in  these  cells  to  such  a  low  level  that  it  impedes  the  ability 
to  remove  water,  which  could  lead  to  the  rapid  decrease  in  cell 
voltage  and  the  eventual  failure  of  the  whole  stack. 

6.  Conclusions 

The  GDM  intrusion  in  PEM  fuel  cells  were  investigated  both 
analytically  and  experimentally.  The  mechanical  behaviors  of 
various  GDM  were  characterized  in  compressive,  flexural,  and 
shear  tests  and  used  in  a  finite  element-based  numerical  model 
to  calculate  the  channel  intrusion  of  GDM.  To  measure  the 
intrusion  ex  situ  and  to  validate  the  model  calculation,  a  test 
fixture  with  a  1  mm  slot  was  devised  to  measure  and  to  simulate 
the  intrusion  of  GDM  in  a  flow  field  with  1  mm- wide  channel. 
Intrusion  of  six  types  of  GDM  from  four  manufacturers  were 
measured  and  compared  with  the  model  calculation.  Excellent 
agreement  was  found  among  these  samples  in  spite  of  a  wide 
range  of  GDM  type  and  mechanical  property,  suggesting  the 
material  tests  adapted  in  this  study  could  be  useful  tools  for 
GDM  manufacturers  in  product  development  and  quality  con¬ 
trol.  To  investigate  the  effect  of  GDM  intrusion  on  the  fuel  cell 
performance,  a  simplified  flow  redistribution  model  was  devel¬ 
oped,  which  showed  that  a  20%  flow  reduction  could  be  induced 
if  reactant  gas  channels  were  to  have  an  additional  5%  intrusion 
beyond  that  of  the  nominal  ones. 

Combining  the  flow  redistribution  model  and  the  GDM  intru¬ 
sion  model,  a  parametric  study  was  conducted  to  investigate  the 
impact  of  GDM  material  properties  and  flow  field  geometry  on 
the  reactant  flow  distribution.  The  study  showed  that  for  flow 
fields  with  a  smaller  channel  depth,  the  reactant  gas  flow  could 
be  very  sensitive  to  the  variation  in  GDM  thickness,  compress¬ 
ibility,  shear  modulus,  and  modulus  of  elasticity.  Within  the 
measured  variation  in  a  commercially  available  GDM,  it  was 
found  that  the  reduction  of  reactant  gas  flow  could  range  from 
10.5%  for  the  flow  field  with  1  mm-deep  channels  to  46%  for 
the  one  with  0.25  mm-deep  channels. 

When  a  set  of  GDM  from  two  different  production  lots  was 
used  in  the  performance  testing  of  a  30-cell  fuel  cell  stack,  it  was 


found  that  the  cells  with  the  GDM  lot  of  higher  intrusion  not 
only  had  lower  cell  voltage  but  also  higher  cell-to-cell  voltage 
variation.  It  was  concluded  that  the  GDM  intrusion  and  intrusion 
variation  had  contributed  to  the  anomalous  behavior  of  this  fuel 
cell  stack.  The  results  from  this  study  showed  that  the  product 
variation  in  GDM  would  need  to  be  more  tightly  controlled  in 
order  to  produce  robust,  efficient,  and  compact  automotive  fuel 
cells. 
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